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Abstract. 

I review the current status of theoretical studies of galaxy formation. I out- 
line the importance of the physics of baryonic component in galaxy formation by 
showing results obtained by using two major tools, semi-analytical approaches 
and cosmological simulations. In particular, I emphasis on roles of feedback in 
galaxy formation and discuss whether apparent confhctions between the stan- 
dard theory of structure formation, the cold dark matter model, and observa- 
tions can be solved by the feedback. I also discuss future prospects in numerical 
simulations of galaxy formation. 



1. Introduction 

Understanding galaxy formation is a challenging problem whose solution will re- 
quire a concerted approach combining observational, semi-analytical and numer- 
ical work. Perhaps the most important progress in this area is the establishment 
of cold dark matter (CDM) model. The CDM model has steadily gained accep 



tance since it was first m ooted in the early 1980s (jPeeblesll 19821 ; iBlumenthal et al 



19841 ; iDavis et al.l Il985l ) . This model has a tremendous predicting power and 



many of these predictions have turned out to be successful. Strong supports 
for this model come from the meas urement of temperature anisotrqpies i n the 
cosmic microwave background (e.g. ISmoot et al.lll992l : ISpergel et al.l 20071 ) and 



two l arge galaxy survey s: the two-degree Field Galaxy Redshift S urvey (2dF- 
GRS IColless et al.ll200lh and the Sloan Digital Sky Survey (SPSS lYork et~al 



2000l ). which have permitted the most ac curate measurements to date of the 



power spectrum of galaxy c l ustering (e .g. iPercival et al.l 120071; iTegmark et al 



20041 . 120061 : iPope et alJ [20041 : ICole et al.1 120051 : iPadmanabhan et al.ll2007l b This 



model specifies the values of the fundamental cosmological parameters and gives 
the initial condition for structure formation of the Universe such as galaxies. We 
are now able to study galaxy formation by ab initio approaches based on the 
CDM model where the structure formation is characterized by the hierarchical 
clustering. 

The second key advance that makes progress in understanding galaxy for- 
mation possible is the unveiling of the high redshift universe. Observations of 
galaxies over a range of redshifts allow us to compare their properties at dif- 
ferent epochs in the history of the U niverse. The unprece dented faint imaging 
of galaxies by the Hubble Deep Field ([Williams et al.l 19961 ) , comb ined with the 



Lym an-break dropout technique to isolate high redshift galaxies (jSteidel et al 



1996), was essential in making possible the first determination of the cosmic 
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star formation history ( Madau et al.|[l996l ; ISteidel et al.l Il999|). The detection 



of the emission from galaxies at sub-m illimeter wavelengths ( Smail et al J 1 19971 : 



Barger et al.lll998l ; iHughes et al.lll99 8) offers the chance of uncovering up heavily 



dust obscured galaxies which are too faint to appear in optical surveys. 

While it is the observations that has led the progress in this field, it is the 
theory that is needed to connect observed snapshots of evolving galaxy popula- 
tion. Our lack of understanding of the key physics underpinning galaxy forma- 
tion, such as star formation and feedback, enforces us to introduce some phe- 
nomenological models. I will illustrate the importance of such baryonic physics 
and how we are beginning to get to grips with it. 



2. Formation of dark matter halos 

Initial tiny density fluctuations grow owing to gravitational instability and even- 
tually form virialized objects, called dark matter halos. In a CDM universe, the 
process of perturbation growth is dissipationless. This means that the total 
kinetic energy of a system of dark matter retained, although energy can be con- 
verted from potential to kinetic. Hence, studying the growth of dark matter 
structure in a CDM universe is essentially straightforward to do by cosmological 
iV-body simulations. 



2.1. The abundance of dark matter halos 

The first at tempt to calculate the ab undance of gravitationally bound objects 
was mad by iPress" k Schechterl ([19741 ) by assuming a Gaussian density field and 
smoothing the field on different scales. The mass function predicted by this 
simple calculation is a grees surprisingly well with the results obta i ned from N- 
body simulations (e.g. lEfstathiou et al.lll988l : lLacev fc Cold [19941 ) . ISheth et"aL 
(120011 ) presented a model i n which they replace d the spherical collapse model 



with an ellipsoidal collapse. iJenkins et al.l ( 200ll ) established the mass function 
of dark matter halos using a suite of iV-body simulations and proposed a fitting 
formula that encapsulates the numerical results . This fitting formul a was found 
to produce good fitting to down to ~ 1O 1O M ( Springel et al.l [20051 ). 

By utilizing the halo mass function, one can derive a mass-to-light ratio that 
guarantees a match between the theoretical prediction and the observed group 
luminosity function. A mass-to-light ratio that guarantees a match between the 
theoretical prediction and the observed group luminosi ty function is a strong 
function of halo mass (lYang et al.ll2003l : lEke et al]l2004T ). It is lowest for halos 
If mass ~ 1O 12 /i _1 M0 and rises by a factor of ~ 6 to lower and higher mass. 
To understand this mass dependence of galaxy formation, we should consider 
baryonic physics. 



2.2. Density profiles of dark halos 

The increasing computer power and the advent of new simulation techniques 
allow to study the formation of individual dark matter halos in full cosmolog- 
ical context (i.e. from CDM initi al conditions). Multi-mass simulation tech- 
niques, so-called zoom simulations. ( Navarro h Benzlll99l1 : iKatz White! 19931: 
Frenk et al.lll996l ) resolve individual dark halos with ~ 10 4 particles ( Navarro et al 
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19961. 11997ft, or even ~ 10 6 particles (|Moore et al.1 Il998l : iKlvpin et afl Il999al : 
Okamoto fc Habdll999ft . 

NFW found that the density profiles of dark halos follow a universal law: 



p(r) 



Po 



(r/r s )(l + r/r s 



(1) 



where r s is a characteristic radius. The profile is characterized by the inner 
slope a = din p / din r = —1 for r < r s and the outer slope a = —3 for r S> r s . 
While many authors have confirmed that the NFW profile provides good fit, 
simul ations with higher resolution form halos with inner slopes as st e ep as a = 
-1.5 (|Moore et al.lll999bl : lFukushige & Makinoll200lh . IHavashi et aD (|2004ft has 
claimed that there is no universal slope for the inner density profile and the inner 
slopes are steeper than NFW and shallower than Moore's, i.e. —1.5 < a < —1 
for r <C r s . 

The existence of such a universal density profile makes strong predictions 
on the shape of the rotation curve of galaxies that can be confronted against 
observations. Such comparisons have revealed that the inner slope of the den- 
sity profiles of a < —1, seems too steep and inconsistent with rotation curves 

studies of dark matter dominated galaxies, i.e. low surface brightness galaxies 

dMcGaugh fc de Blokl[l998l : iFlores k Primackill994l : iMoore et al.lil999bl : Ide Blok k Bosma 
20021 ). Some authors have argued that the rotation curves are consistent with 
CDM models if one considers the triaxiality of dar k matter halos (jHavashi k Navarro 
2006; IValenzuela et al.ll2007l : IHavashi et aDl2007l ). However we should consider 
effects of baryonic component properly since the presence of a galaxy disk at the 
center of a dark matter halo makes dark halo significa ntly rounder than dark 
matter halos formed in dark matter only simulations ( Bailin et a l. 2005) and 
adiabatic contraction due to disk formation makes dark matter halos more con- 
centrate as we shall see later. The problem is not only at the very center of low 
surface bright galaxies but a l so seen at intermedia te radii of any disk galaxies 
([Navarro k SteinmetdbOQObl ; iMcGaugh et al.ll2007l ). This problem in the CDM 
density profile is called the cuspy core problem. 



2.3. Substructure within dark matter halos 



Until the late 1990s, cosmological iV-b qdy simulations suggested that dark mat - 
ter halos were smooth and featureless ( Summers et al. 19951 ; Frenk et al.lll996l ). 
Higher resolution simulations which typically contained more than 10 6 particles 
in each dark matter halo showed that this phenomenon, called over-merging, was 
a numerical artifact (iGhigna et al.|[l99^ ; IMoore et al.lll998l ; IKlvpin et al.lll999a1 ; 
Okamoto fc Habdll999ft . Once a halo enters within the virial radius of a more 
massive halo, it is referred to as a satellit e halo or substructure within the larger 
halo {sub-halo). lOkamoto k Habel ( 19991 . [2000l ) provided a method to construct 
merging hi stories of individual sub-h alos in order t o follow their formation and 
evolutions. IMoore et al.l ( 1999al ) and IKlvpin et al.l ( 1999b! ) pointed out that the 
number of satellite halos in a Milky Way-sized halo is an order of magnitude 
larger than the number of the satellite galaxies of the Milky Way. This problem 
is called the satellite problem. 
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3. Semi-analytic models 



A powerful way of studyin g hierarchical galaxy formation is semi-analytic mod- 
elling of galaxy formation. White k, Reesl ( T978I ) proposed that galaxy formation 
was a two stage process, with dark halos forming in dissipationless, gravitational 
collapse, with galaxies forming inside these structures, following the radiative 
cooling of baryons. They also argued the importance of an additional process, 
fe edback, for avoiding the production of more faint galaxies than are observed. 



White Sz Frenki (j 199 II ) produced a galaxy formation model that included many 



of today's models: C DM, radiative coo l ing of gas, star formati o n, fee dback, and 
stellar populations. iKauffmann et al.1 ( 19931 ) and ICole et ail ( 19941 ) presented 
the first models to track the formation and evolution of galaxies in the set- 
ting of evolving dark matter halos. Below I briefly describe procedures that 
are commonly used in semi-analytic models. A more c omprehensive review on 
semi-analytic galaxy formation models can be found in lBaughl (|2006l l. 



3.1. Basic ingredients 

Dark halo merger trees Semi-analytic models require the formation history of 
each dark matter halo, commonly called the merger tree, in which a galaxy 
forms and evolves. Merger trees can be constructed using a Monte-Carlo ap- 
proach by sampling the distrib ution of progenitor masses predicted by the ex- 
tended Press-Sc hechter theory ( Lacev fc Cold [19941 ; ISomerville fc Kolattl [l999l ; 
ICole et al.ll2000h . More faithful way to construct merger trees is extracting from 
iV-body simulations which have sufficiently frequent outputs. T he first attempt 
to ext r act merger trees from iV -body simulations was made by iRoukema et al 



( 19971 ). IKauffmann et all ( 19991 ) s tudied galaxy distribution usin g merger trees 
generated by iV-body simulations. lOkamoto &: Nagashima (|200ll ) utilizes merger 
trees of sub-halos obt ained in high-resolution cosm ological TV-body simulations 
of clusters of galaxies ( Okamoto Habdll999l . 20001 ) to investigate the morphology- 
density relation of cluster galaxies ( Dressier 198Cll ). Th is powerful method was 
firstly used for semi-analytic stu dies of cluster galaxies ( Okamoto Nagashima! 
120011 . 120031 : ISpringel et atll200lh and has become a standard technique for semi- 



analy tic models combined with iV-body simulations ( Croton et al.|[2006l : Bower et al 



analy 
2006) 



Gas cooling The basic model of how ga s cools inside dark matter halos was 
set out in detail bv IWhite fc Frenki (|l99ll ). Semi-analytic models assume that 
when a dark matter halo forms, gas contained in the halo is shock heated to the 
virial temperature of the halo: 



1 fArrip 

2 k 



V c , 



(2) 



where fi is the mean molecular weight of the gas, m p is the mass of a hydrogen 

atom, k is Boltzmann's constant, and V c = (GM v ; r / V v i r )2 is the circular velocity 
of the halo with mass M v i r and virial radius r v ; r . This gas with T = T Y - ir is 
called the hot gas. It is also assumed that the hot gas distributes parallel to 
the dark matter. The relation between the circular velocity and halo mass is a 
function of redshift and cosmology. Gas can subsequently cool from the hot halo. 
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The cooling rate of the gas is dependent on the temperature and metallicity of 
the gas. Hence a cooling time can be specified by dividing the thermal energy 
density of the gas by the cooling rate per unit volume: 



WW= = /(A(p g ,T vir ,Z)), (3) 




where p g is the gas density which is a function of radius and the function A is 
a cooling function. The cooling radius, r coo i is computed as the radius where 
the cooling time is equal to the time-step used in the merger tree. The gas 
within r coo i accretes to the center of the halo at rate set by the cooling time (or 
dynamical time if the dynamical time is longer than the cooling time) keeping 
its angular momentum. The cooled gas is assumed to have ~ 10 4 K at which 
the atomic cooling becomes inefficient and this gas is called the cold gas. 

Star formation The lack of a theory of star formation enforces and allows a 
simple estimate of the global rate of star formation in a model galaxy: 

M* = (4) 

T* 

where the star formation rate, M*, depends on the amount of cold gas available, 
M co i(j and a star formation timescale, r*. The star formation timescale is usually 
represented as 



(5) 



where r° is either the dynamical time T(j yn = r ga i/t! ga i within the galaxy, which 
is a function of redshift, or some constant and the dependence on the circular 
velocity is introduced with two free paramteters, K and a*, to r eproduce the ob- 
serve d cold gas fractions in spirals as a function of luminosity ( Cole et al J 1 19941 . 



2000). The star formation timescale that does not depend on the dynamical 
timescal e (i.e. redshift) has been supp orted from analysis on the formation of 



qua sars 



3 (i.e. reasniit) nas been supported ironi analysis on tne lormation oi 
( Kauffmann &: Haehneltj |2000| ) . the evolution of damped Lyq system s 



( Somerviile et al.ll2001 ). and the galaxy number counts ( Nagashima et al. 200ll ). 



Since the value of a* is usually assumed to be less than -1, the star formation 
timescale is longer for smaller galaxies. This naturally explains the galaxy down- 
sizing by keeping the gas in small galaxies until today, though the physical origin 
of the form of the star formation timescale is unclear. 

Merging of galaxies When two or more dark halos merge to form a new halo, 
the hot gas components in the progenitor halos are immediately merged and 
constitute a hot gas component in the new halo. In contrast, galaxies do not 
merge immediately. The central galaxy of the main progenitor halo is defined 
as the central galaxy of the new halo and the rest of galaxies are regarded 
as satellite galaxies. Satellite galaxies sink to the halo center owing to the 
dynamical friction and merge to t he central galaxy. The satell ite galaxies could 
merge each by random collision ( Somerviile Primack 1999). If the merging 



two galaxies have similar mass, the merger is considered as the major merger 
and a starburst occurs. The new galaxy becomes a pure bulge galaxy. In the case 
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of a minor merger, the disk of the larger galaxy is not destroyed and a merging 
satellite is absorbed into either the disk or the bulge component of the larger 
galaxy. Note that so me authors introduced small sta rbursts induced by minor 
mergers. For example lOkamoto fe Nagashima (2003) claimed that a starburst 



induced by a minor merger is necessary component in order to form the galaxy 
population with intermediate morphology. 

Feedback Feedback processes arguably have the largest impact on the form of 
the theoretical predictions for galaxy properties, whilst at the same time being 
amongst the most difficult and controversial phenomena to model. The most 
common form of feedback used in models is the ejection of cold gas from a 



•li 

galac tic disk by a supernova (SN) driven wind (e.g. iLarsonl 119741 ; iDekel Silk 



1986). The reheated cold gas could be blown out to the hot gas halo, from which 



it may subsequently recool (sometimes called retention feedback) , or it may even 
be ejected from the halo (ejection feedback) and left unable to cool until it is 
incorporated into a more massive halo later. The ejection rate of the cold gas 
by SN driven wind is parameterized as follows: 

/ y \ ahot ■ 

where Vhot an d aohotp are free parameters. Considering the fact that the poten- 
tial wells are shallower in smaller galaxies, the value of cthot should be negative. 
From energy balance, Qhot = —2 is obtained. However many authors used 
much stronger dependence on the circular velocity (i.e. aw < —2) in order 



to reproduce the faint end slope of luminosity functions (e.g. iCole et al.l 11994 
Nagashima fc Yoshiill200ll ) 



iBenson et al.l (|2003l ) carried out a systematic study of the impact of various 
feedback mechanisms on the form of the predictions for the galaxy luminosity 
function. The standard SN driven winds, while helping to reduce the number of 
faint galaxies to match observations, were found to overproduce bright galaxies. 
Stronger feedback of this mode would tend to weaken the break in the predicted 
luminosity function rather than enhance it by wiping out galaxies around L*. 
The thermal conduction and the ejection mode of feedback (superwinds) did 
the better job. However they had to assume an unphysically high conductivity 
in the conduction model and an i mplau sibly high efficiency of fe edback in the 
superwind model. ICroton et al.l ( 20061 ) and iBower et al.l ( 20061 ) implemented 



simple AGN feedback schemes into their models, respectively. By assuming the 
radio mode feedback from AGN quenches gas cooling in quasi-hydrostatic hot 
halos, they reproduced the break in the luminosity function. 

Chemical evolution Chemical evolution of gas and stars in a galaxy is impor- 
tant because (i) th e cooling rate of the gas, A (T,Z), strongly depends on the 
gas metallicity (e.g. ISutherland Dopita1ll993l ); (ii) the metallicity with which 



stars are born has an impact on the luminosity and colors of the stellar popula- 
tion; and (iii) the optical depth of a galaxy, which determines the extinction of 
starlight due to dust, scales with the mtallicity of its cold gas. As stars evolve, 
they return material to the interstellar medium (ISM) with an enhanced metal- 
licity as the form of stellar winds or SN explosions. The amount of gas returned 
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to the ISM per unit mass of stars forms is there for dependent on the form of 
the initial mass function (IMF) of the stellar population. 

The first attempt to follow the chemical evolution of galaxies in semi- 
analytic models considered type II SNe, using instantaneous recycling approxi- 
mation. While the yield of metals is determined by the chosen IMF, the effective 
yield of metals depends on a number of differing metallicities such as (i) the mix- 
ing of hot gas reservoirs of differing metallicities by halo mergers, (ii) accretion 
of gas by cooling to the disc cold gas which could have different metallicity, and 
(iii) feedback processes which deliver metals in the cold gas to the hot halo gas. 

Recently type la SNe have been included in semi-analytic models. iThomas 

(1999) was the first to consider the delayed enrichme nt due to SNe la by usin g 
star formation histories extracted from the model of iKauffmann et al.l (I1993T). 
but n eglecting any inflow or outflow of gas and metals. iNagashima fc Okamotd 
(2006) produced the first semi-analytic model which self-consistently integrated 
the impact of SN e la, using the simplificat ion of assuming a fixed time delay for 
SN la explosions. INagashima et ail (2005a, b) performed the first fully consistent 
calculation including both type II and la Sne in the semi- analytic models. 



3.2. Successes and failures 



Semi-analytic models now successfully re produce the presen t day galaxy lumi- 
nosity functions by invoking spurewind s dBenson et al. 20031 ) or AGN feedback 
(ICroton et al.ll2006l : iBower et al]l200d ). H ower et al.l also reproduced luminos- 
ity functions at higher redshifts. However readers should bear in mind that any 
phenomenon which leads to a model successfully matching the observed break 
in the luminosity function uses up a dangerously high fraction of the energy re- 
leased by supernova explosions. We also have to wait to see whether the models 
which invoke AGN feedback can account for the observed AGN activities. 

Semi-analytic models have a difficulty to match the zero-point of the Tully- 
Fisher relation and the normalization of the luminosity function at the same 
time. The difficulty lies on the fact that the CDM halos are too centrally con- 
centrated and therefore it might be a genuine problem of the CDM. Scaling 
relations such as the Tully-Fisher relation rely on the calculation of the scale 
lengths of galaxies. Semi- analytic models assume that hot halo gas has the 
same specific angular momentum as its host halo and the angular momentum 
is conserved when gas becomes cold and accretes to a disk. However, numerical 
simulations have shown th at angular mom e ntum of gas can be transferred to 
the dark matter. Moreover, Okamoto et al. ( 20051 ) showed that the direction of 
the angular momentum vector of the main progenitor halo changed significantly 
duri ng galaxy formation . 



Baugh et al.l (|2005l ) were able to reproduce the number counts of sub- 



millimeter selected galaxies and the luminosity function of Lyman-break galax- 
ies, whiles retaining a fair match to the present day optical and far infrared 
luminosity functions. However, this was only possible with the controversial 
assumption of a top-heavy IMF for star formation in merger driven starbursts. 
This model can al so account for the metal licity of th e intracluster med i um and 
elliptical galaxies ( Nagashima et al.ll2005al lbl). though INagashima et al.l ( 2005b! ) 
failed to reproduce the observed trend of a/Fe increasing with the velocity 
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dispersion of the galaxy. So far there is no explanation for this trend in the 
framework of the hierarchical galaxy formation. 



4. Gasdynamic simulations 



Another powerful tool to study galaxy formation in cosmological context is gas- 
dynamic simulations. In the gasdynamic simulations gas dynamics is solved 
based on governing equations. This is the major advantage of the gasdynamic 
simulations over the semi-analytic models in which gas distribution is assumed. 
There are two principle algorithm in common use to follow the hydrodynamics 
of gas in an expanding universe: particle-based Lagrangian schemes, which em- 
ploy a technique called smoothed particle hydrodynamics (SPH) (e.g. Monaghanl 



grid-based Eulerian schemes (e.g. iRvu et al 



1992 : ICouchman et~aDl995l : [P earce et al.lll999|;ISpringel &: Hernquisti 



19931 : ICen fc Ostrikerl 



2003ft . 



and 
19991). To 



date, most of gasdynamic simulations of galaxy formation have been peformed 
by SPH, because its Lagrangian nature is sutable to follow the collapse of gas 
clouds in cosmological volume. 

Prescriptions used in gasdynamic simulations in order to treat baryonic 
physics such as star formation and feedback are surprisingly similar to those 
used in semi-analytic models. The main difference is that such processes are 
formulated based on the local quantities in simulations rather than global quan- 
tities such as the circular velocity of the halo. For example, the rate of change 
of the internal specific energy of i-th particle by cooling is calculated as 



A(pi,Ti,Zi) 



Pi 



(7) 



where Uj, pi, T i} and Z{ are the specific internal energy, density, temperature, 
and metallicity of the z-th particle, respectively. Note that in Eq.Q, the virial 
temperature of the halo and the mean metallicity of the hot halo gas are used. 
Although m any simulations followed the chemical evoluti on, only a handful of 
studies (e.g. iKawata &: Gibsonll2005l ; lOkamoto et al.ll2005l ) employed the metal- 
licity dep endent cooling function and other use the cooling rate of the primordia l 
gas (e.g. lAbadi et al.l 120031 : ISommer-Larsen et al.l 120031 : iRobertson et al.l 120041 ; 
iGovernato et al. 20071 ) in spite of the fact that the cooling rate of the gas with 
the solar metallicity is more than an order of magnitude higher than that of the 
primordial gas. Usually, the inverse Compton cooling and effects of time evolving 
uniform ultraviolet background are also considered in the energy equation. 

Dense (p > p s f ~ 0.1 cm~ 3 ) and cold (T ~ 10 4 K) gas is eligible to form 
stars. Star formation rate of the i-th particle is given by 



C* - 



oc p 



1.5 



(8) 



where C* is the dimensionless s tar formation eff iciency parameter. To match 
the observed Kennicutt rela tion ( Kennicu tt 1998), C* ~ 0.05 is often used (e.g. 
Navarro fe Steinmet d EoOOal ) . 

Each star particle represents a single stellar population. Thus evolved stars 
explode as SNe. Instantaneous recycling is often assumed. Some authors relaxed 
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Figure 1. Stellar distributions at z = 0. The no-burst, density-burst and 
shock-burst models are shown from left to right. The edge-on and face-on 
views of stars are given in upper and lower panels, respectively. The brightness 
indicates the projected mass density, and the same scaling is used for each 
model. All of these galaxies are obtained from the same initial conditions. 



this approximation and cal culated the number of SNe and returned mass as a 



function of the SSP's age (|Okamoto et al.1 1200.1 120071 : iGovernato et alJl2007l ). 



They also included SNe la. The feedback energy is distributed to the surrounding 
gas. It was found that depositing energy in the form of thermal energy had 
almost no e ffect because the energy was quickly radiated away in dense regions 
(|Katdll992l ). Note that since the cooling function depends on the temperature, 
not only the amount of energy given to the gas but also the total mass of gas 
to which the energy is given affect the thermal evolution. For example, giving 
all the feedback energy to the nearest one particle has much stronger effect 
than giving it to 50 particles. To circumvent this problem, feedback energy is 
often given by the form of the kinetic energy (e.g. iNavarro Sz Steinmet3 2000a: 
Springel HernquistJl2003l )or cooling for the particles that receive the feed back 
energy is shutting-off for a while (~ 10 Myr) ( Thacker Sz Couchman 200 ll ). 
A self-consistent AGN feedback model w as firstly incorporated into cos- 



mological simulations of galaxy formation by ( Okamoto et al. 20071 ) in which 



they discriminated the quasar mode a nd radio mod e base d on the mass accre- 
tion rates onto the central blackholes. Siiacki et al. ( 20071 ) showed that such a 
discrimination successfully produced bright red galaxies. 



4.1. Formation of spiral galaxies and effects of feedback 

The first attempts to simulate the formation of a spiral galaxy from CDM initial 
conditions generally failed, producing objects with overly centrally concentrated 
gas and stars. It became immediately apparent that the root cause of this prob- 
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lem was a net transfer of angular momentum from the baryons to the dark mat- 
ter halo during the aggregation of the galaxy thro ugh mergers ( Navarro Sz Benzl 
19911 : iNavarro fe Whitdll994l : lNavarro et alJll995l ). This is known as the angular 
momentum problem. It was suspected from the start that its solution was likely 
to involve feedback processes that would regulate the supply of gas to the galaxy. 

More recent simulati ons within the CDM f r amew ork have produced more 
promising disk galaxies. iThacker Couchman ( 200ll ) obtained a disk galaxy 
having reasonable size by assuming that gas heated by SN explosions can adia- 
batically expand without rad iative loss of energy fo r a wh ile. There simulation, 
however, stopped at z = 0.5. ISteinmetz fc Navarrol ( 20021 ) did not assume such 
strong feedback and found that a broad range of galaxy morph ologies could 
be pro duced by gas accretion and galaxy mergers. In related work, lAbadi et al 



( 200 J) obtained a disk galax y which resembled observed early-type disk galaxies 
bommer- Larsen et all (120031 ) were also able to generate a variety of morphologi- 
cal types, including disks, by assuming high star formation effi ciency and strong 
feedba ck at high redshift to prevent early collapse of baryons. Governato et al.l 
( 20071 ) used the feedback model similar to IThacker Sz Couchman ( 200ll ) in or- 
der to match the Tully-Fisher relation. There Tully-Fisher relation was, how- 
ever, calculated from rotation speed of stars younger than 4 Gyr-old, which 
is much older populations than those observed by Ha. In fact, the rotation 
speed of the cold gas i n their galaxy was too fast compared with observations. 
Robertson et al.l ( 2004! ) adopted a multiphase model for the ISM which stabi- 
lized gaseous disks against the Toomre instability, and produced a galaxy with 
an exponential surface brightness profile but insufficient angular momentum. 
These simulations, however, did not take the metallicity effects on cooling into 

acco unt. 

I Okamoto et al. ( 20051 ) also assumed a multiphase model for the ISM and a 
top-heavy IMF for stars formed in starbursts, as required for semi- analytic mod- 



els to match the number counts of bright submilimeter galaxies (jBaugh et al 



2005) and metal conten ts in the intracluster medium and elliptical galaxies 
([Nagashima et alj|2005al[bf). By varying the criteria for starbursts in their sim- 
ulations, lOkamoto et"ail l2005) ware able to produce galaxies with a variety of 
morphological types, from ellipticals to spirals, starting from exactly the same 
initial conditions (Fig. [I]) . 



The density profiles In Fig. [21 density profiles of the host dark matter halos of 
the three galaxies in Fig. Q] are shown. The profile obtained by the simulation 
without baryonic component (Ob = 0) is also shown as the DM-only simulation. 
It is clear that the dissipatinal nature of the baryons makes the inner density 
profile steeper. The no-burst galaxy which has the weakest feedback effects has 
the steepest profile. 



The evolution of angul ar momentum Zavala et al. ( 20071 ) analyzed two of the 



galaxies simulated by lOkamoto et al.l (|2005l ). the no-burst galaxy (left panel 
in Fig. [T]) and the shock-burst galaxy (righ). The no-burst galaxy is referred 
as the bulge-dominated galaxy and the shock-burst galaxy is called the disk- 
dominated galaxy. They traced the time evolution of the angular momentum 
of three Lagrangian components: the dark matter particles lie within the virial 
radius at z = (halo), the 10% most bound dark matter particles within the halo 
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Figure 2. The left panel shows density profiles of dark halos for the DM 
only, no-burst, density-burst, and shock-burst simulations by the solid, dotted, 
dashed, and dot-dashed lines, respectively. The right panel is the same as the 
left one but for the total mass density profiles. 



(inner halo), cold gas and stars within 0.1r v i r (galaxy). The tidal torque theory 
predicts that the amplitude of the angular momentum of a Lagrangian region 
evolves in proportion to a 3 / 2 until the region reaches the maximum expansion. 
Then the region starts to collapse and is decoupled fr om the tidal field, and thus 
the angular momentum stays constant (|Whitdll984l ). 

The angular momentum evolution of the halo component (upper-left panel 
of Fig. [3]) clearly shows two distinct phases, indicated by the dashed (L oc a 3 / 2 ) 
and dotted (L = const.) lines. The inner halo shows quite different evolution 
(upper- left panel). Since this subregion has a higher overdensity than the halo as 
a whole, it reaches maximum expansion earlier. Until z = 0, the inner halo loses 
90% of the angular momentum by mergers. Lower panels in Fig. [3] show angular 
momentum evolutions of galaxies (cold baryons) for the bulge-dominated galaxy 
(left) and disk-dominated galaxy (right). The evolution of the bulge-dominated 
galaxy (lower-left panel) is surprisingly similar to that of the inner halo (upper- 
right). This implies that baryons that are contained in the bulge-dominated 
galaxy are concentrated in the dense subsystems of dark matter at high redshift 
and lose their angular momentum as the inner halo loses its angular momentum 
during the collapse. On the other hand the evolution of the disk-dominated 
galaxy (lower-right) resembles that of the halo as a whole (upper-left). Strong 
feedback from the stellar population with a top-heavy IMF avoids early collapse 
of baryons and therefore baryons are diffusely distributed as the hot halo gas. 
When this gas cools, it accretes to the disk keeping its angular momentum. 
Thus, the suppression of the early collapse of baryons by some form of feedback 
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Figure 3. Upper- left: the evolution of the rms radius (top panel) and the 
specific angular momentum (middle and bottom panels) of the dark mater 
particles that lie within r vlr at z = 0. In the bottom panel the dashed line 
indicated L oc a 3 / 2 and the dotted line L = const. Upper-right: same as 
the upper-left but for the particles that makes the inner dark matter halo. 
Lower- left: Evolution of the rms radius (upper panel) and specific angu- 
lar momentum evolution (lower-panel) for the baryonic component of the 
bulge-dominated galaxy. Lower-right: same as the lower-left but for the disk- 
dominated galaxies. All physical quantities are normalized to their values at 
the present day. 



is a key to solve the angular momentum problem. Inter estingly, the model use d 
in the disk-dominated galaxy (shock-burst simulation in lOkamoto et all ( 20051 )) 
was also able to repro duce the luminosity function of the Milky Way satellites 



(|Libeskind et al.ll2007l ) 



5. Conclusions 

I have presented results on recent efforts to model the formation and evolution 
of galaxies in hierarchically clustering universes using two complementary tech- 
niques: semi-analytic models and gasdynamic simulations. Feedback processes 
are the most important ingredient in these models which have reproduced the 
observed luminosity function and solved the angular momentum problem and 
the satellite problems. It is interesting that the feedback recipes that can solve 
the angular momentum problem also solve the satellite problem simultaneously 
(jLibeskind et al.1 120071 : iGovernato et al.1 120071 ). However in both techniques a 



dangerously high fraction of the energy released by SNe is used up. Sometimes 
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more than 100% of energy is used in simulations. An excuse for doing so is 
that the ISM has a multiphase structure and therefore the cooling should not 
be as effective as in the simulations where the multiphase structure is not re- 
solved. This uncertainty can be easily removed if we use the resolution that is 
sufficiently high to resolve star forming gas with density as high as bh ~ 100 
cm 3 and temperature as low as ~ 100 K. Simulations of the ISM on the galactic 
scales have recentl y started fTask er fc Brvanl l2006; Ro bertson Kravtsovll2007l ; 



Saitoh et al.1 12007m . ISaitoh et alJ "if2007) found that in these high resolution sim- 



ulations star formation is not controlled by the star formation time-scale defined 
in equation ([8]) but regulated by the time scale of gas supply form low density 
regions to high density regions, which is p/p ~ 5 x tdyn- While there is still long 
way to go, resolving multiphase structure of the ISM in cosmological simulation 
will be a crucial step to make our physically motivated models to models based 
on physics. The radiat ion feedback has be en largely ignored or only crudely 
included in our models. iMurrav et al.l ( 20051 ) suggested that large scale galactic 



winds seen around local and high redshift starburst galaxies are driven by radi- 
ation pressure rather than energy from SN explosions. To test this hypothesis, 
we have to carry out radiation hydrodynamics in galaxy formation simulations. 
Semi-analytic models also need refinements. While they have already have many 
parameters and are well-complicated, they are still too simple to take dynami- 
cal effects such as angular momentum transfer into account. Close interactions 
between simulations and semi-analytic models will be required for further im- 
provements and better understanding of galaxy formation. 
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